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The algebraic method for computing the complete point symmetry group of a system of differential 
equations is extended to finding the complete equivalence group of a class of such systems. The 
extended method uses the knowledge of the corresponding equivalence algebra. Two versions of 
the method are presented, where the first involves the automorphism group of this algebra and the 
second is based on a list of its megaideals. We illustrate the megaideal-based version of the method 
with the computation of the complete equivalence group of a class of nonlinear wave equations with 
applications in nonlinear elasticity. 


1 Introduction 

Classes of (systems of) differential equations are (systems of) differential equations that contain 
parameters, which may be constants or functions, and are collectively referred to as arbitrary 
elements. Most systems of differential equations of science and engineering involve such arbitrary 
elements and hence are classes rather than single systems. The wide occurrence of classes of 
differential equations is the reason why their study has a prominent place in group analysis of 
differential equations m [261 EH ED]. The classification of Lie symmetry properties of systems 
from a class depending on the values of the arbitrary elements is called the group classification 
problem for this class. 

Equivalence transformations of a class of differential equations are point (resp. contact) trans¬ 
formations in the underlying space of independent and dependent variables, involved derivatives 
and arbitrary elements that map each system from the class to another system from the same 
class [261 EQj. Such transformations play a central role in group classification of differential 
equations since they can be used to select the simplest representative among similar systems. 
Continuous equivalence transformations form a connected Lie group, which can thus be found 
using the infinitesimal method by computing the corresponding Lie algebra, the so-called equiva¬ 
lence algebra, i.e., the set of vector fields which generate one-parameter subgroups of equivalence 
transformations. Due to its algorithmic nature, this is a standard way for computing equivalence 
transformations in papers on group classification of differential equations mm- 

The main shortcoming of the infinitesimal method is that it misses discrete equivalence trans¬ 
formations. That is, not the complete equivalence group is obtained but only its continuous com¬ 
ponent. The situation is completely analogous to the construction of the maximal Lie invariance 
group for a single system of differential equations, which is carried out within the framework 
of the infinitesimal method; but in order to present the corresponding complete group of point 
symmetries, also discrete symmetry transformations have to be computed. 

Finding complete point symmetry groups (and, analogously, complete equivalence groups) is 
a nontrivial task if the direct method is applied In the direct method one 

computes point symmetries without the detour of infinitesimal generators. The direct method 
hence also allows one to find discrete point symmetries and thus to construct the complete point 
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symmetry groupie This method does not use the linearization of the determining equations, 
which enables the algorithmic computation of Lie symmetries, and thus usually one has to solve 
a system of nonlinear partial differential equations; which is a challenging endeavor. For systems 
with a nontrivial continuous component of the complete point symmetry group, a rehnement 
for the pure direct method can be used, the so-called algebraic method. The algebraic method 
is simpler in that it allows one to restrict the principal form of point transformations before 
invoking the direct method. In particular, the solution of nonlinear systems of partial differential 
equations can often be completely avoided. 

The central idea of the algebraic method for the computation of discrete symmetries of 
differential equations is the following: Given a system of differential equations C, let qc denote 
its maximal Lie invariance algebra, which is known. Push-forwards of vector fields defined on the 
corresponding space of independent and dependent variables by point symmetries of C induce 
automorphisms of Qc- This property can be used in order to restrict the form of those point 
transformations that can be symmetries of C. The restricted form can then be substituted into 
the system £ from which point one still proceeds with the direct method. In fact, the algebraic 
method just simplifies the application of the direct method. 

The original version of the algebraic method was proposed by Hydon nanaiis]; see also [TO] . 
It is based on the immediate use of the above property of point symmetries of C and re¬ 
quires the explicit computation of the entire automorphism group of Qc] this is why we call 
it the automorphism-based version of the algebraic method or, shortly, the automorphism-based 
method. The computation of the entire automorphism group is a comprehensive or even im¬ 
possible task for infinite-dimensional Lie algebras. Therefore, for a proper application of the 
automorphism-based method the algebra Qc has to be hnite dimensional, which is a severe 
limitation of the applicability of this version of the algebraic method. 

There are countless examples of systems of differential equations with inhnite-dimensional 
maximal Lie invariance algebras, especially in hydrodynamics. This is why in [3] an essential 
refinement of the algebraic method was proposed. It is built around the notion of fully char¬ 
acteristic ideals (shortly called megaideals); these are vector subspaces of a Lie algebra g that 
are invariant under any transformation from the automorphism group of g mElEH], see also 
[121 Exercise 14.1.1]. Usually, the requirement for a linear nondegenerate mapping of g to leave 
invariant each of the megaideals of g produces enough constraints to understand the rough struc¬ 
ture of automorphisms. Moreover, in [11IIII29] several propositions were proved that enable the 
systematic computation of new megaideals using known ones. Given a system of differential 
equations £, let Qc be its maximal Lie invariance algebra. If a sufficient number of megaideals 
of the maximal Lie invariance algebra Qc of the system £ is known, the condition of their invari¬ 
ance with respect to automorphisms of Qc can be used in order to derive maximal restrictions 
on the form of a general point symmetry transformation of £ that are obtainable within the 
framework of the algebraic method. This is the essence of the megaideal-based version of the 
algebraic method or, shortly, the megaideal-based method. 

The effectiveness of this version has been demonstrated by computing the complete point 
symmetry group of the barotropic vorticity equation on the beta-plane [3], the vorticity equa¬ 
tion on the sphere [7], the quasi-geostrophic two-layer model |3| and the system of primitive 
equations [S]. For each of these (systems of) differential equations, the corresponding maximal 
Lie invariance algebra is infinite dimensional. 

The same techniques as applied to compute the complete point symmetry group of a system 
of differential equations can also be used to compute the complete equivalence group of a class of 
systems of differential equations. Moreover, the megaideal-based method is even more effective 

^Discrete symmetries like alternating signs of independent and/or dependent variables can easily be found ad 
hoc, by empiric methods. Even if such trivial symmetries exhaust a complete set of discrete symmetries that are 
independent up to composing with continuous symmetries, which is often the case, the problem is to prove that 
this is really the case. Such a proof cannot be realized within an empirical approach. 
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here than for symmetry groups as equivalence groups are often infinite dimensional. Showing 
how to realize the computation of equivalence groups using the algebraic method is the purpose 
of the present paper. Both versions of the algebraic method are discussed. 

We illustrate the megaideal-based method by the computation of the equivalence group for 
the class of nonlinear wave equations of the general form 


Utt = f{x, Ux)Uxx + g{x, Ux), 

(la) 

with the nonvanishing condition 


ifux^ 9uxUx) 7^ (0)0)- 

(lb) 


This condition is essential as it guarantees that all equations from the class ([T]) are really non¬ 
linear, cf. [2]. Physically, this class arises, in particular, in the study of nonlinear elasticity. 

The study of the class ([T]) plays a prominent role in the field of group classification of differen¬ 
tial equations. It was initially considered in where one-dimensional symmetry extensions of 
the kernel algebra with respect to a fixed finite-dimensional subalgebra of the infinite-dimensional 
equivalence algebra of this class were found. In other words, a partial preliminary group classi¬ 
fication problem was solved. A similar direction was taken in m , where the partial preliminary 
group classification of the class ([T|) with respect to one-dimensional subalgebras of an infinite¬ 
dimensional subalgebra of the equivalence algebra was considered. The group classification 
problem for systems of the form vt = a{x, v)wx and wt = b{x, v)vx related to equations from 
the original class ([T]) was first investigated in [TU] and then exhaustively solved using the alge¬ 
braic method of group classification in [20] . A study related to the group classification problem 
of the class ([T]) was carried out in [IB], resulting in the construction of second-order differen¬ 
tial invariants of the equivalence algebra. Finally, the complete group classification problem 
for the class © was solved in |2] by the enhanced algebraic method of group classification, 
thus completing the over twenty years long investigation of the point symmetry properties of 
this class. 

The further structure of the present paper is as follows. In Section [21 we recall some defini¬ 
tions related to classes of systems of differential equations, their equivalence algebras and their 
equivalence groups. Section [3| contains the necessary material explaining the algebraic method 
for the computation of (complete point) equivalence groups. Both versions of the method— 
automorphism- and megaideal-based—are introduced and presented in the form of step-by-step 
procedures. In Section [H we illustrate the megaideal-based method with the computation of 
the equivalence group of the class ([T]) . In the final Section O we summarize our results and give 
some thoughts about future directions of investigation. 


2 The equivalence group and equivalence algebra 

In this section we recall the central definition of the equivalence group and the equivalence 
algebra of a class of (systems of) differential equations. These are two of the fundamental 
concepts that play a key role in group classification of differential equations. As such, they are 
discussed extensively in the papers and books on this subject, see, e.g., [DElElESlEZlEHllso]. 
It is natural to first rigorously define the notion of a class of differential equations. 

Let Cq: L(x, 0(q)(a:,U( p))) = 0 denote a system of differential equations parameterized 
by the tuple of arbitrary elements 0(x,U(p)) = (0^(x, M(p)),..., 0^(x, U(p))). Here and in the 
following, the tuple x = (xi,..., Xn) consists of the independent variables and includes both 
the tuple of dependent variables u = {u^, ..., u™) as well as all the derivatives of u with respect 
to X of order up to p. By 0(g) we denote the partial derivatives of the arbitrary elements 6 of 
order not exceeding q for which both x and U(p) act as the independent variables. 
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In the definition of the class m we explicitly included the inequality / ( 0 , 0 ), 

which restricts the values that the two arbitrary elements / and g of the class ([T]) can take. Such 
a set of so-called auxiliary conditions on the arbitrary elements 0 is a necessary and important 
component of the precise definition of the class of differential equations but it is often posed in 
an implicit way. The procedure can be formalized as follows. 

Consider the system of auxiliary differential equations S{x, 9(^qi){x, U(^p))) = 0 and inequal¬ 
ities S(x, / oH in which again both x and act as independent variables. 

Components of the tuples S and S are smooth functions of x, U(^p) and 9(qi)- It is thus required 
that the tuple of arbitrary elements 9 runs through the solution set, denoted by 5, of both the 
auxiliary equations S = 0 and inequalities S / 0. 

The above discussion can be summarized in the following definition: 

Definition 1. The set {Cq | 0 E 5} denoted by is called a class of differential equations 
defined by parameterized systems Cq and the set S of arbitrary elements 9. 

In order to define the equivalence group of the class TI 5 it is convenient to first introduce 
the notion of admissible transformations and the notion of equivalence groupoid [21 [271 [281130]. 
Given the arbitrary elements 9,9 £ S with the associated systems Ce and Cg from the class C\s, 
the set T(0, 9) of point transformations that map the system Cq to the system Cg is called the set 
of admissible (point) transformations from Lq to Cg. In particular, the maximal point symmetry 
(pseudo)group Gq of the system Cq coincides, by definition, with the set T{9,9). 

A triplet {9,9,(j)), where 9,9 £ S with T{9,9) 7 ^ 0 and f £ T(^)i^)) is called an admissible 
transformation in the class In other words, an admissible transformation is a triplet con¬ 
sisting of a pair of similar systems (a source and a target system) and a point transformation 
connecting these two systems. 

The partial binary operation of composition is naturally defined for pairs of admissible trans¬ 
formations for which the target system of the first admissible transformation and the source 
equation of the second admissible transformation coincides, {9,9,(j))o{9,9,(j)) = (9,9, (pof). Each 
admissible transformation is invertible, {9,9,cj))~^ = {9,9,(f>~^), i.e., the inversion of admissible 
transformations is a unitary operation defined everywhere. The set of admissible transforma¬ 
tions of the class £[5 that is endowed with the operations of composition and taking the inverse 
is called the equivalence groupoid of this class and denoted by = Q""{C\s)- All the groupoid 
axioms are obviously satisfied for . 

For a general class of differential equations, it is not possible to relate its equivalence groupoid 
with a group of point transformations that act in the extended space of (x,M(p), 0 ), respect the 
contact structure on the space of (x,M(p)) and preserve the class. However, it makes sense to 
single out a maximal part of the equivalence groupoid that admits such a relation; then the 
associated group is referred to as the equivalence group. 

Definition 2. The (usual) equivalence group, denoted by = ^^(£ 15 ), of the class £[5 is 
the (pseudo)group of point transformations in the space of {x,U(^p-^,9), each element of which 
satisfies the following properties: It is projectable to the space of (x,M(p/)) for any p' with 
0 ^ p' ^ p. The projection is the p'th order prolongation of <h|(a;,u). For any 9 from S 

its image ^9 also belongs to S. Finally, £ T{9,^9). 

Here we say that a point transformation ip: z = ip{z) in the space of variables z = (zi,..., Zk) 
is projectable on the space of variables z' = (zq,..., Ziy) with 1 ^ ii < • • • < ^ /c if the 

expressions for the transformed variables z' depend only on z!. The projection of ip to the 
£-space is denoted by ip\zi, z' = ip\zi{z'). 

Elements of are called equivalence transformations of the class C\s- Each equivalence 
transformation induces a family of admissible transformations parameterized by the arbitrary 

^Other kinds of inequalities (>, <, etc.) as well as collections of them are also relevant. 
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elements, {(0, <I>0, | 6 G 5}. If the entire equivalence groupoid is induced by the 

equivalence group in the above way, then transformational properties of the class are 
particularly nice and the class is called normalized 121 EH |28l [30]. The notion of normalized 
classes of differential equations serves as basis for developing the algebraic method of group 
classification of differential equations. 

If the arbitrary elements 9 do not depend on derivatives of u of order greater than p', p' < p, 
then the group can be assumed to act in the space of (x, tt(p'), 9). 

In some cases, e.g., if the arbitrary elements 9 depend on x and u only, we can neglect 
the condition that the transformation components for (x, u) of equivalence transformations do 
not involve 9, which gives the generalized equivalence group of the class £[5 

[lllETllso]. Each element d> of jj is a point transformation in the (x, u, 0)-space such that 
for any 9 from S its image ^9 also belongs to S, and <b(-, •, 9{-, •))|(a;^„) € T(0, ^9). 

It was mentioned in the introduction that the computation of point transformations of sys¬ 
tems of differential equations using the direct method usually involves the solution of nonlinear 
systems of differential equations. Thus, rather than studying finite equivalence transformations, 
it is common in the framework of group analysis of differential equations to only consider their 
infinitesimal counterparts. This restriction linearizes the nonlinear determining equations for 
finding such transformations, making their computation essentially algorithmic. Thus, to the 
continuous component of the equivalence group one may associate the Lie algebra of 
vector fields in the space of (x,M(p), 0 ), which for any 0 ^ p' ^ p are projectable to the space 
of {x,U(^pi-j) with the property that their projections to the space of {x,U(^pi)) are the p'th order 
prolongations of their projections to the space of {x,u). Moreover, these vector fields are the 
generators of one-parameter groups of equivalence transformations of the class £[ 5 . The Lie 
algebra is called the (usual) equivalence algebra of the class £[ 5 . 

In a similar way, by considering the generalized equivalence group of the class £[5 instead 
of the usual one and by allowing the components of vector fields that correspond to x and u to 
also depend on 0 , one can define the generalized equivalence algebra of the class £[ 5 . 

The (usual or generalized) equivalence algebra g~ of a class £[5 can be found by invoking an 
infinitesimal invariance criterion completely analogous to that used for the computation of the 
infinitesimal generators of one-parameter point symmetry groups of a single system of differential 

equations mm- 

In order to use the infinitesimal invariance criterion properly, one has to consider a vector 
field Q on the extended space of independent variables x, derivatives and arbitrary ele¬ 
ments 9. The components of Q corresponding to derivatives of u with respect to x should 
respect the contact structure of the space of independent and dependent variables. In other 
words, these coefficients should be obtained using the general prolongation formula [25l l26] . 
The vector field Q should be prolonged to the derivatives of the arbitrary elements assuming 
both X and as the independent variables; again, the general prolongation formula should be 
used for this purpose. The prolonged vector field is then applied to the joint system A of the 
general form L(x, = 0 of systems from the class £[5 and of the auxiliary conditions 

defining the set <S of values of the arbitrary elements 9. If necessary, differential consequences 
of equations from the system A have to be considered. The remaining computational procedure 
then parallels the computation of infinitesimal generators of Lie point symmetries and involves 
the solution of an overdetermined linear system of partial differential equations. 

3 Algebraic method 

As was mentioned in the introduction, any point symmetry transformation T of a system of 
differential equations £ generates an automorphism of the maximal Lie invariance algebra of £ 
via push-forwarding of vector fields in the space of system variables. This condition implies 
constraints for the transformation T which are then taken into account in further calculations 
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using the direct method [a 0 0 na la la ES]. The set of transformations found in the 
way described constitute the complete point symmetry group of the system C including both 
continuous and discrete point transformations. 

The above algebraic method can be easily extended to the framework of equivalence trans¬ 
formations. A basis for this is given by the following simple theorem, which is proved in the 
same way as the similar assertion for point symmetries. 

Theorem 1. Let C\s he a class of (systems of) differential equations, and the equiva¬ 
lence group and the equivalence algebra of this class (of the same type, namely, either usual or 
generalized onesjH Any transformation T from induces an automorphism of gT via push¬ 
forwarding of vector fields in the relevant space of independent variables, derivatives of unknown 
functions and arbitrary elements of the class. 

Proof. Consider an arbitrary vector field Q G g~. The local one-parameter transformation group 
Gq = {exp(e(5)} associated with Q is contained in G"^. Then the one-parameter transformation 
group Gq = {Texp(eQ)T~^}, which is similar to Gq with respect to T and whose generator is 
T^tQ, is also contained in This means that the vector field T*Q belongs to An arbitrary 
push-forward respects the Lie bracket of vector fields, [T*Q, %Q'] = %[Q, Q'] for any Q, Q' € g~. 
Therefore, 7^ is an automorphisms of g~. □ 


In other words, for each element T of there exists an automorphism T of g"^ such that 
the condition T*Q = TQ is satisfied for any Q from g~. 

It should be noted that while the correspondence T ^ 71 defines a representation of 
in Aut(g~), this representation is often unfaithful. This is always the case if the equivalence 
group (resp. the equivalence algebra g~) has a nontrivial center. Moreover, the representation 
image Gf is a subgroup of the automorphism group Aut(g~) that may be essentially smaller than 
the entire automorphism group Aut(g~). In other words, there may exist many automorphisms 
of g~ that are not induced by elements of 

The discussion on continuous point symmetries of a single system of differential equations C 
given in [7] can be adopted to equivalence transformations as well. In particular, continuous 
equivalence transformations of the class C\s can be composed and thus generate a connected 
normal subgroup, denoted by ?7~, of the equivalence group G^. The transformations from 
induce mappings on the equivalence algebra g"^ that are internal automorphisms of g~ thus 
generating the normal subgroup Int(g'^) of the automorphism group Aut(g~). So, Int(g~) = Uf. 
The representatives of the elements of the factor group G'" in will be interpreted as 
discrete equivalence transformations of the class C\s. It is clear from the construction that 
it makes sense to consider only discrete equivalence transformations that are independent up 
to composing them with continuous equivalence transformations of the class C\s. When 
is infinite dimensional, the factorization is usually difficult but it is sometimes still possible to 
realize it and thus to identify the discrete equivalence transformations. An example will be given 
in Section 01 

Remark 1. Classes of systems of differential equations with a single dependent variable (i.e., 
with m = 1) may possess contact equivalence transformations that are not prolongations of 
point transformations. In this case the entire theory can be extended to contact equivalence 
transformations. 

®The further consideration does not depend on what kind of equivalence groups and equivalence algebras (usual 
or generalized) are involved. For consistency, we should just handle the usual (resp. generalized) equivalence 
algebra when looking for the usual (resp. generalized) equivalence group and impose relevant a priori restrictions 
on the admitted structure of equivalence transformations. 
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3.1 Version based on automorphisms 

Theorem [T] allows us to extend the automorphism-based method for finding point symmetry 
groups of systems of differential equations [laiiaiis] to equivalence groups of classes of such 
systems. The extension properly works if the corresponding equivalence algebra 0 "^ is of finite 
nonzero (and, moreover, low) dimension since the knowledge of the entire group Aut( 0 ~) is 
needed. A new feature is that one should consider appropriate point transformations and vector 
fields in the relevant extended space of independent variables, derivatives of unknown functions 
and arbitrary elements, i.e., of (x,n(p),0). If the arbitrary elements 9 do not depend on deriva¬ 
tives of u of order greater than p', p' < p, then these appropriate objects can be restricted to 
the space of (x, 9). In the case p' > 0 they should be consistent with the contact structure 

of the space of In the course of computing usual equivalence groups, they should also 

be projectable to the space of the system variables {x,u). 

We present the extension of the automorphism-based method in the form of a step-by-step 
procedure. In what follows the indices i, j and k run from 1 to n = dim 0 ~, and we assume 
summation over repeated indices. 

1. Let T |5 be a class of (systems of) differential equations, find the equivalence algebra 0 ~ 
of this class by the infinitesimal invariance criterion. Suppose that n = dim^^ < oo. 

2 . Fix a basis {Qi ,..., Qn} of 0 ~, and compute the structure constants of 0 ~ in this basis, 

[Qi,Qj] = d^jQk- 

3. The general form (ap of automorphism matrices in the basis {Qi ,..., Qn} is obtained via 

solving the system of algebraic equations c^ijia\ 0 

4. Take the general form of an appropriate point transformation T in the space of x, rt(p) and 
9 and push-forward the vector fields Qi,, Qn- Then we set %Qi = alQj, i = 1,... ,n, 
where (aj) is the general form of automorphism matrices. Consequently equating the 
corresponding vector-field components of right- and left-hand sides of these equations 
produces a system of differential equations for components of the transformation rl 

5. Integrate the system and use the obtained intermediate form of T within the framework 
of the direct method in order to complete the system of constraints for T and to produce 
the final form for T. 

Recall once more that the difference in the application of the above procedure to the compu¬ 
tation of the usual (resp. generalized) equivalence group is that is the usual (resp. generalized) 
equivalence algebra; additionally, for the usual equivalence group one restricts oneself to trans¬ 
formations that are projectable to the space of x and u while for the generalized equivalence 
group one takes T to be a general point transformation in the space of (x, u, 9) without any 
additional a priori restrictions. 

Remark 2. The procedure described above can be modified by factoring out the group of 
continuous equivalence transformations of the class which can easily be computed once the 

Automorphism groups have been computed for many finite-dimensional Lie algebras, in particular for all 
semi-simple Lie algebras m and Lie algebras of dimension not greater than six, see [S] (9] 1101129| and references 
therein. The automorphism group of a decomposable Lie algebra is easily constructed from the automorphism 
groups of the decomposition components [9]. Finding the automorphism group of a Lie algebra g can be simplified 
if the chosen basis of 0 is consistent with a known megaideal hierarchy of g. 

®The system derived involves parameters of the automorphism group Autfg^). Therefore, in the course of 
integrating the system one in fact needs to solve a compatibility problem in order to find values of the parameters 
for which the system is consistent with respect to components of T. If the group G7 does not coincide with 
Autfg^), which is often the case, then there exist values of the parameters for which the system is inconsistent. 
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equivalence algebra g~ is known. As ?7~ = lnt( 0 '^), in order to do the modification, instead 
of considering the general form of automorphism matrices in Step 3, one can factorize general 
automorphisms by inner automorphisms. More precisely, suppose that for each equivalence class 
of automorphisms with respect to lnt( 0 '^) one can take a representative in such a way that these 
representatives constitute a matrix group. Then, this group is isomorphic to the factor group 
Aut( 0 ~)/Int( 0 ~). In general, the form of representatives is more restrictive than the general 
form of automorphism matrices, which may simplify the equations for T derived in Step 3 of the 
above procedure. The realization of the modified procedure results in the discrete equivalence 
transformations of the class C\s composed with elements of the center of ?7~. 

For an equivalence algebra of large finite dimension or, especially, for an infinite-dimensional 
equivalence algebra the computation of its automorphism groups may be a difficult problem, 
which also needs additional efforts for its rigorous formulation in the case of infinite dimension. 
Moreover, a significant part of automorphisms may not be realized via push-forwarding of vector 
fields by appropriate point transformations in the relevant space of independent and dependent 
variables, involved derivatives and arbitrary elements. In such a case, the exhaustive description 
of automorphisms is not too principal as then the equivalence transformations induce only 
a (small) proper subgroup of the automorphism group. 

3.2 Version based on megaideals 

It has been pointed out in Section [TT] that for the automorphism-based method, it is crucial that 
the equivalence algebra is low dimensional. If this is not the case, the computation of the 
automorphism group Aut( 0 '^) becomes an intricate problem. It is a quite common situation for a 
class of differential equations to be studied that its equivalence algebra g~ is infinite dimensional 
and thus a more appropriate version of the algebraic method is needed. Mathematically, the 
main problem with the version reviewed in Section 13.11 is that in computing the automorphism 
group Aut( 0 ~) one also has to find the elements Aut( 0 ~) \ which cannot be used for the 
computation of 

This deficiency is remedied by working with megaideals of 0 ~ rather than with Aut( 0 ~) itself 
HEIIH]. That is, we employ the fact that G~i C i for any megaideal i of 0 ~, rather than invoking 
the condition C Aut( 0 ~)|^ 

We recall that a fully characteristic ideal (or, shortly, megaideal) i of a Lie algebra g is a vector 
subspace of g that is invariant under any transformation from the automorphism group Aut( 0 ) 
of 0 [11[29], cf. [T2l Exercise 14.1.1]. We thus have that Ti = i for each megaideal i of g, whenever 
T is a transformation from Aut( 0 )Ji| Megaideals of 0 are ideals and, moreover, characteristic 
ideals of g. In the present context, the Lie algebra in this definition is the equivalence algebra g”^ 
of the class of differential equations C\s- 

In order to make the version of the megaideal-based method most effective, it is essential to be 
able to construct wide sets of megaideals. The more megaideals are known, the better restrictions 
on the form of admitted automorphisms can be derived. While it might be a complicated problem 
to obtain a complete list of megaideals of g”^, for the practical application of the megaideal-based 
method the more tractable problem of finding a set of megaideals that gives maximal restrictions 
on point transformations should be tackled. Thus, it is for example not essential to consider 
megaideals that are sums of other megaideals, since they give weaker constraints than their 
individual summands jointly. 

^Although the stronger condition GJ'i = i is always satisfied, the inverse inclusion G~i D i is trivial in view of 
the presence of the identical automorphism in G^, gives no constraints on elements of G~ and, therefore, can be 
neglected. 

^The invariance condition C i of i with respect to implies i C Ti, which gives, together with the 

condition Ti C i, the equality Ti = i. 



In [H[7l|8l|29] it was shown how new megaideals can be computed from known ones. For the 
sake of completeness of the present paper, we collect the statements on megaideals from these 
papers. The central observation for practical computations is that many megaideals of a Lie 
algebra 0 can be constructed without prior knowledge of the automorphism group Aut( 0 ) using 
the following obvious assertions: 

Proposition 1. Let ii and i 2 be megaideals of a Lie algebra 0 . Then we have: 

1. The improper subalgebras of q (i.e., the zero subspace and 0 itself) are megaideals of q. 

2. The sum ii + i 2 , the intersection ii n i 2 and the Lie product [ii,i 2 ] of megaideals are 
megaideals. 

3. If i 2 is a megaideal of ii and ii is a megaideal of 0 then 12 is a megaideal of 0 . Thus, 
megaideals of megaideals are again megaideals. 

4 . The centralizer (resp. the normalizer) of a megaideal is a megaideal. 

5. All elements of the derived, upper and lower central series of a Lie algebra are its megaide¬ 
als. It thus follows that the center and the derivative of q are megaideals. 

6. The radical r and nil-radical n (i.e., the maximal solvable and nilpotent ideals, respectively) 
of g as well as different Lie products, sums and intersections involving 0 , r and n f] 0 ,tr], 
[r, r], [ 0 ,n], [r, n], [n, n], etc.) are megaideals of g. 

In order to find more megaideals without computing automorphisms, we also apply an as¬ 
sertion, which in general has no clear interpretation in terms of distinguished object related to 
the structure of the corresponding Lie algebra. It was first proved in [7], but for the sake of 
reference we repeat the proof here. 

Proposition 2. If io, ii and i 2 are megaideals of 0 then the set s of elements from io whose 
commutators with arbitrary elements from ii belong to i 2 is also a megaideal of g. 

Proof. It is clear that the set 5 is a linear subspace of 0 . Consider an arbitrary element zq 5. 
Hence zq G io and [zq, zi] € 12 for any zi G ii. Then for any T G Aut( 0 ) and any zi G ii we have 
Tzq € io and \fZzQ,zi] = [TzO)TT“^zi] = ‘I[zq,‘I~^zi] G i 2 as‘I~^zi G ii and thus [zo,T~^zi] G i 2 . 
This means that Tzq G s. Therefore, s is a megaideal of 0 . □ 

Assertion (iv) of Proposition [1] is a particular case of Proposition [2l where io = 0 and i 2 = {0} 
for the centralizer of ii, and io = 0 and i 2 = ii for the normalizer of ii. 

Remark 3. If m is a finite-dimensional megaideal and its automorphism group Aut(m) is already 
known, all megaideals of m can be found by direct computation according to the definition of 
megaideals as subspaces of m invariant with respect to Aut(m). In the course of calculating the 
automorphisms we can use knowledge about structural megaideals of m such as the center, the 
radical, the nilradical, elements of the derived, lower central and upper central series of m and 
different megaideals related to structural megaideals via certain operations. 

Theorem [1] implies that any transformation T from satisfies the condition 7j=i = i for 
each megaideal i of 0 ~. In the course of the derivation of constraints for T, this condition is 
interpreted in the following way: Let the megaideal i consist of the vector fields where 7 
runs through a parameter set T,, and be spanned by the vector fields Qy parameterized by 7 ' 
running through a subset T- of Pi. The condition Tli = i implies that for any 7 ' G T', there 
exists a 7 G Ti such that %Qy = Q-y. 

Similarly to the automorphism-based method, the megaideal-based counterpart can also 
be split in a few algorithmic steps. Some steps are the same or almost the same as for the 
automorphism-based method. Even so, for the convenience of further application we describe 
each step in its entirety. We also re-interpret the condition for megaideal invariance in a manner 
more convenient for practical use. 
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1 . Given a class of (systems of) differential equations C\s-, find the equivalence algebra 

of this class by the infinitesimal invariance criterion. Let = {Q.^ I 7 G r} for some 
parameter set L. 

2. Fix a set {Qy \ 7 ' € F'} of vector fields spanning g~ = {Q^i \ 7 ' G F') with F' C F, 
and compute commutation relations for all pairs of spanning vector fields. 

3. Using Propositions [Hand [21 Remark[3|and other tools, construct as wide a list of megaideals 
of the equivalence algebra g~ as possible. Optimize the list by the exclusion of inessential 
megaideals, which give no new constraints for automorphisms of g~ as compared with other 
megaideals. In particular, megaideals being sums of other megaideals are not essential. 

4. Take the general form of an appropriate point transformation T in the space of x, U(^p'j 
and 6 and push-forward the vector fields Qy, G Fh For each 7 ' we choose, from the 
above list, the minimal megaideal \^i = {Q^ | 7 G F^, C F} containing Qy and set 
71) Qy = Q-yy, where the parameter satisfies the constraints singling out the set Fy, 
from the entire set F. Consequently equating the corresponding vector-field components 
of right- and left-hand sides in the equation %Qy = for each 7 ' leads to a system 
of differential equations for components of the transformation rH The solution of this 
system results in an intermediate form for T. 

5. Use the derived form of T within the framework of the direct method in order to complete 
the system of constraints for T and to produce the final form for T by solving the completed 
system. 


4 Complete equivalence group of the class 
of nonlinear wave equations 

In order to clarify all the algorithmic steps, in this section we demonstrate in detail the compu¬ 
tation of the usual equivalence group of the class of nonlinear wave equations ([U using the 
megaideal-based method. The preliminary presentation of this result is contained in the arXiv 
preprint of [2]. 

4.1 Equivalence algebra 

The equivalence algebra of class © including both linear and nonlinear equations was first 
computed in HZ]. It was then shown in |2], that the class consisting of only nonlinear equations 
admits the same equivalence algebra. Here we use both the representation of the equivalence 
algebra and the notation introduced in [2]. The equivalence algebra g~ of the class © is spanned 
by the vector fields 

= udu + Uxdu, + gdg, = tdt - 2fdf - 2gdg, = dt, 

— ‘fdx '-Px'^xdux. T 2,ipxfdf 4“ y^xxUx f dg, (2) 

Gi-ip) = i>du + 'Gxdu, - il^xxfdg, = tdu, = t^du + 2dg, 

where (p = (p{x) and ijj = 'ip{x) run through the set of smooth functions of x. It has been pointed 
out in [2| that while the elements of the equivalence algebra of class © are vector fields in the 

^Similarly to the automorphism-based method, the system obtained involves the parameters 7.^/, 7' G F'. 
Therefore, integrating the system requires solving a compatibility problem in order to find values of the parameters 
for which the system is consistent with respect to components of T. In addition to that the group G7 may be 
a proper subgroup of Aut(0~), for the megaideal-based version there may be one more reason for the system to 
be inconsistent for some parameter values: In general, even the complete megaideal hierarchy cannot properly 
capture delicate constraints for elements of Autfg'^). 
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space of (i, a:, ti( 2 ),/, 5 ), for practical purposes it suffices to present only their projection to the 
space of {t, x,u,Ux, f, g)- The reason for this is that the vector field components corresponding 
to ut, Ux, utt, utx and Uxx-, can be obtained via prolongation from the components corresponding 
to t, X and u. At the same time, we explicitly include the components associated with the 
derivative Ux in the representation of vector fields ([ 2 ]) spanning since the arbitrary elements / 
and g depend on Ux- Therefore, for the proper computation of commutation relations, the 
components with du^ are crucially needed. 

For the construction of megaideals, it is convenient to also recall the nonvanishing commu¬ 
tation relations between the vector fields of g~. They read 

[gw,v'^] = GW, 

[V\T^] = 

[p\V^]=V\ [V\T^]=G{1), = 

[v{ip),g{ip)] = g{ipii)x)- 

Therefore, steps 1 and 2 of the megaideal-based method were in fact realized in ElE]. 

4.2 Megaideals of equivalence algebra 

In order to compute the complete equivalence group of the class ([I|) by the megaideal-based 
method (cf. Section I4.3jl , we need to know a set of megaideals of the equivalence algebra of 
this class. 

Let g = g”^ for the sake of notational simplicity. Using Proposition [U it is easy to compute 
the following megaideals of g~: 

q' = {V\V{ip),g{^l^),F\T^), q" = {V{ip),g{^l^),F^), g'" = (PM,e(V’)), 

C0(g'") = {V\V\g{l),F\F^), C,,{q"') = {V\g{l),F\F^), 

c,,{Q") = {g{i),F\F^), z,„ = {g{i),F^), Zg, = (^(1)), 

Rg = {V^,V\v\g{il^),F\F^), Rg.. = (aW), 

where a', Rq, Zq and Ca(b) denote the derivative, the radical and the center of a Lie algebra a 
and the centralizer of a subalgebra b in a, respectively. We present proofs only for the last two 
equalities. 

The linear span Si = {V^,g{'4>),F^ ^F"^) obviously is a solvable ideal of g. Moreover, 
it is the maximal solvable ideal of g. Indeed, suppose that Si C i and i is an ideal of g. Then 
there exists a smooth function ^ of x which does not identically vanish such that the vector field 
P((j) belongs to i. As i is an ideal of g, for an arbitrary smooth function tp oi x the commutator 
[D{C),'D{ip)] = V{Cipx — Cxf) belongs to i. If ( is not a constant function, we define the following 
series of operators; 

j,k = 1,2,.... 

It is possible to prove by induction that / 0, j, /c = 1, 2,... . Moreover, 

as G i, we have € i^-^\ j,k = 1 , 2 ,... , i.e., i^-^^ 7 ^ { 0 } for any nonnegative integer j. 
This means that the ideal i is not solvable. If C is a constant function, we can set (" = 1. We 
choose any smooth function (p of x with ipxx ^ 0 and denote ipx by C- As the commutator 
[D{l),'D{ipy\ = T>{0 belongs to i, the consideration for the previous case again implies that the 
ideal i is not solvable. Therefore, Si is really the maximal solvable ideal of g, i.e., Rg = Si. 

The linear spans 2 = (^(V')) is an Abelian and, therefore, solvable ideal of q'”. The maximality 
of this solvable ideal is proved in the same way as for Si. Hence Rg/// = S 2 . 

The same megaideals can be obtained in different ways. For example, (I?(l)) = Zg/ = Zg///. 
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To find one more megaideal which will be used in the course of the computation of the 
complete equivalence group of the class ([I]) by the algebraic method in Section 14.31 we should 
apply a more sophisticated technique than before. There are two ways to do this. 

The way directly based on the definition of megaideals is to calculate the automorphism 
group Aut(m) of the finite-dimensional megaideal m = Cg(g"') = (P*, 7 ^*, and then 

determine megaideals of m as subspaces of m which are invariant with respect to Aut(m). In the 
course of calculating the automorphisms we use the knowledge about simple megaideals of m, 

The presence of the above set of nested megaideals is equivalent to that for any automorphism A 
of m, its matrix in the basis {g{l),T ^is upper triangular with nonzero 

diagonal elements. In particular, 

A'P^ = 014^(1) + a24-T^ + 034 + 0447 ^*, 

AV ^ = 015^(1) + a - 2 b ^^ + 035-^^ + <245'^* + a55^^ 

where 044055 / 0. As \P^,V^] = and A G Aut(m), we should have [AV^,AV^] = AV^. 
Collecting coefficients of basis elements in the last equality, we derive a system of equations with 
respect to o’s which implies, in view of the condition 044 7^ 0, that 055 = 1, 034 = 0, 024 = 044035 
and oi4 = 044025 — 045024. As 034 = 0, we get that the span {V^ is a megaideal of m 
and, therefore, of 0~. 

The other way for finding the megaideal ( 7 ^^, 17 ( 1 ), which is based on Proposition [21 
allows us to avoid the above computation of automorphisms of the megaideal m and, thus, 
is much simpler; Choosing io = ii = Cg/(0'") = ( 7 ^*, ^( 1 ), and i2 = Zg/ = ( 17 ( 1 )) in 

Proposition O we obtain the megaideal s = {V^,g{l),F^). 


4.3 Calculation of equivalence group 


The careful study of megaideals of the equivalence algebra 0 ~ in Section 14.21 according to the 
third step of the megaideal-based procedure supplies us with a sufficient store of megaideals in 
order to commence the computations directly concerned with equivalence transformations. 

Theorem 2. The equivalence group C of the class © consists of the transformations 

C2Ux + 


t = Clt -|- Co, X = ¥^(x), U = C2U + C4,l + c^t -|- V'(x), 


Ux = 




/ = %/- ~9=^ 


, C2Ux +'lpx r I f , n 
C 2 g H- TxxJ - iPxxJ + 2C4 

Tx 


( 3 ) 


where cq, ..., C4 are arbitrary constants satisfying the condition C1C2 / 0 and g? and fj run 
through the set of smooth functions of x, ipx 7 ^ 0 . 


Proof. The group consists of nondegenerate point transformations in the joint space of 
variables t, x and u, the first derivatives ut and Ux and the arbitrary elements / and g^ which 
are projectable to the variable space and whose components for first derivatives are defined via 
the first prolongation of their projections to the variable space. Thus, the general form of a 
transformation T from is 


t = T{t,x,u), x = X(t,x,u), u = U{t,x,u), 

Ul = U\t, X, U, Ut, Ux), Ux = IX^it, X, U, Ut,Ux), 

f = F{t, X, u, Ut, Ux, f,g), g = G{t, x, u, ut,Ux, f, g). 
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where and are determined via T, X and U and the nondegeneracy condition should be 
additionally satisfied. To obtain the constrained form of T, we will act by the push-forward Tl 
induced by T on the vector fields ([2]) additionally including the terms with du^ and use megaideals 
of the equivalence algebra of the class ([T]) and restrictions on automorphisms of found 
in Section [4.21 The principal part of the consideration is based on objects and properties of 
that are related to the finite-dimensional megaideal m = Cg(g"') = {T>^, ^(1), Recall 

that megaideals being sums of other megaideals are not essential for the computation since they 
give weaker constraints than their summands. For example, the megaideal g" is the sum of g'" 
and Zg/ and hence we do not use it in the further consideration. A list of essential megaideals is 
in fact exhausted by the spans 

(Gil)), {G{l),X^), {G{1),X\J=-^), {V\G{l),T^), {V\V\G{1),X\X^), 

(ew), 

We apply the condition of preserving megaideals following their ordering by inclusion. In other 
words, for each vector field Q from ([2]) we take the megaideal that contains Q and is minimal 
among the above listed megaideals. As a result, we derive an optimized set of constraints for % 
as automorphism of g, 


% gii ) = 011^(1), 

(4a) 

= 012^(1) + 022.F^, 

(4b) 

%F^ = 013^(1) -1- 023^^ -1- azzF'^, 

(4c) 

= 044^(1) -|- 024^^ -|- , 

(4d) 

TFD ^ = 015^(1) -I- 025^^ -I- OssF^ -I- 045^* -|- F*, 

(4e) 

% g {'^) = G { F ^), 

(4f) 

TMF ) = G { r ) + T ){ pn , 

(4g) 


where a’s are constants, 011022033044 7 ^ 0 and ip and (p are arbitrary smooth functions of x. 
The constants o’s completed with 055 = 1 and o^ = 0, 1 ^ i j ^ 5, form a matrix of an 
automorphism of the megaideal m. Tildes over vector fields on the right hand sides of the above 
equations mean that these vector fields are written in terms of the transformed variables. The 
parameter-functions and are smooth functions of x associated with the parameter- 

functions Ip or if, which is indicated by the corresponding superscripts. We will derive constraints 
for T, consequently equating the corresponding vector-field components of right- and left-hand 
sides in each of the conditions @ and taking into account constraints obtained in previous 
steps. As the components of vector fields and of the transformation T associated with the 
derivatives ut and Ux are defined via first-order prolongation involving the similar values related 
to the variables t, x and u, the coefficients of dut and du^ give no essentially new equations in 
comparison with the coefficients of dt, dx and du- This is why we will not equate the coefficients 
of dut and du^. To have a proper representation of the final result, we will re-denote certain 
values in an appropriate way. 

Thus, the equation (Ha|l implies that Tu = Xu = 0, Uu = C 2 and Fu = Gu = 0, where 
the nonvanishing constant an is re-denoted by C2. Then we derive from Eq. (| 4 bp that tUu = 
a22T + ai2, i.e., T = ciI + cq where ci = 02/022 / 0 and cq = -012/022, and Fut = Gut = 0- The 
consequence t^Uu = + 0,2‘iT+ai'i of Eq. (| 4 c|l gives only relations between o’s. In particular, 

“33 = C2 /cp Then the other consequences of Eq. (| 4 cl) are = 0 and Gg = C2 /cf. The essential 
consequences of Eq. () 4 dll are exhausted by Xt = 0, Ut = a24T -|-ai4 and Ft = Gt = 0. Therefore, 
X = p{x) and U = C2 + C4t‘^ + c^t -|- ip^x), where px 7^ 0, C4 = 02404/2 and 03 = 044 -|- 02409. 

As we have already derived the precise expressions for the components of T corresponding 
to the variables (cf. Eq. ([3|)), at this point we could terminate the computation of equivalence 
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transformations by the algebraic method and calculate the expressions for F and G by the direct 
method. At the same time, all the determining equations for transformations from the equiv¬ 
alence group of the class ([T|) follow from restrictions for automorphisms of the equivalence 
algebra 0 ~. This is not a common situation when the algebraic method is applied. Usually it 
gives only a part of the determining equations simplifying the subsequent application of the direct 
method. See, e.g., the computations of the complete point symmetry groups of the barotropic 
vorticity equation and quasi-geostrophic two-layer model in [U Section 3] and [3l Section 4], 
respectively. This is why we complete the consideration of the equivalence group within the 
framework of the algebraic method. 

From Eq. (Hel) we obtain in particular that tUt = a 35 T^-|-a 25 r-|-ai 5 , fFf = F and fGj+gGg = 
G — 035 , where 035 = 2 c 4 /c^ in view of the first of these consequences. 

The equation ([4f| implies the equations 

= ^.Fu^=0, ^^Gu^ - i^a^a^fGg = i^i^F (5) 

The first and second equations of ([5]) are equivalent to = C 2 'ip and F^^ = 0. Then we can 
express the derivative ip'f- via derivatives of {(fx^’xx — fxxi’x), substitute the 

expression into the third equation of ([5]) and split with respect to the derivatives ij^x and ijjxx, 
as the function 'll) is arbitrary. As a result, we obtain F = c^'^ipx'^f and Gu^ = C2(Px^(pxxF, i.e., 
Guj. = C 2 Ci‘^(px^fxxf- The expression for F coincides with the transformation component for / 
presented in the theorem. 

The last essential equation for G is given by Eq. ([3^. Collecting coefficients of dx, du and 
dg in Eq. ([4g]), we have that = ipx^, = ijjx^ and 

^Gx ^xUxGu^ T ‘^0xfGj^ ^xx^x f G g ip--UxF '^XX^I (®) 

respectively. We proceed in a way analogous to the previous step. Namely, we express the 
derivatives (p^- and via derivatives of (p, substitute the expressions into Eq. Q and split 
with respect to derivatives of (p because the function ip is arbitrary. Equating the coefficients 
of (px leads to the equation 

fGf — ipx; (c2Ux T 'Px'pxx 'Pxx^x^' 

The simultaneous integration of all the equations obtained for G precisely results in the 
transformation component for g from the theorem, which completes step 4 of the megaideal- 
based method. 

In order to complete the proof, we should realize the last procedure step, which is in fact 
trivial for the class m- We should just check by the direct computation of expressions for 
transformed derivatives that any transformation of the form ([3]) maps any equation from the 
class ([ll) to an equation from the same class. □ 

After comparing the equivalence algebra and the equivalence group the following 
corollary is evident: 

Corollary 1. A complete list of discrete equivalence transformations of the class 0 that are 
independent up to combining with each other and with continuous equivalence transformations 
of this class is exhausted by the transformations 

{t, X, u, Ux, f, g) {-t, X, u, Ux, /, g), 

{t, X, u, Ux, f, g) ^ ft, -X, u, -Ux, f, g), 

{t, X, u, Ux, f, g) {t, X, -u, -Ux, f, -g). 
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Theorem [2] implies that any transformation T from of the class ([T]) can be represented 
as the composition 

r = 2)'(ci)T*(co)DMD“(c2)T^(c4)3-2(c3)g(V^), 


where the above parameterized equivalence transformations are 


T*(co): 

i 

= t + Co, 

X 

= X, 

u 

= u, 


— '^X'! 

f 

D‘(ci): 

i 

= Clt, 

X 

= X, 

u 

= u, 


— '^X'! 

f 


i 

= t, 

X 

= 

u 

= u, 

^x 


f 

D“(c2 ): 

i 

= t, 

X 

= X, 

u 

= C 2 U, 

^x 

- ^ 2 '^X'! 

f 


i 

= t, 

X 

= X, 

u 

= u + c^t, 

^x 

— "^Xl 

f 

3''(c4): 

i 

= t, 

X 

= X, 

u 

= u + , 

^x 

— '^Xl 

f 

SW: 

i 

= t, 

X 

= X, 

u 

= u + ljj, 

^x 

— tix H“ ' 4 ^X 1 

f 


f, 9 = g, 

-2 r ~ -2 

Cl J 1 9 — Cl g, 

vlf, 9 = 9 + VxxUxf/^Px, 

f, 9 = C2g, 

f, 9 = 9, 

f, 9 = 9 + ‘2 c4, 

f, 9 = 9- i^xxf, 


and the nondegeneracy requires that ciC 2 V 2 a; 7^ 0. These transformations are shifts and scalings 
in t, arbitrary transformations in x, scalings of u, gauging transformations of u with square 
polynomials in t and arbitrary functions of x. 


5 Conclusion 

In this paper we have extended the algebraic method for the computation of complete point 
symmetry groups of single systems of differential equations to the computation of equivalence 
groups of classes of such systems. In general, the equivalence group may include both discrete and 
continuous equivalence transformations. Unlike the infinitesimal method, which aims at finding 
the corresponding equivalence algebras from which the continuous parts of equivalence groups are 
computed, the algebraic method allows also constructing discrete equivalence transformations. 
We have developed two versions of this method, automorphism-based and megaideal-based, 
which may be combined depending on the problem under consideration [7]. The effectiveness of 
both versions rests on the fact that each equivalence transformation induces, via push-forwarding 
the vector helds that constitute the corresponding equivalence algebra an automorphism of 
the algebra g~. This imposes substantial restrictions on the functional form of equivalence 
transformations and strongly simplifies, or even makes trivial, the further application of the 
direct method. 

Each of the versions has own advantages and disadvantages. We briefly recall some of them. 

The automorphism-based version ensures the maximal use of properties of equivalence alge¬ 
bras within the framework of the algebraic method. One more benefit of this version is that 
most continuous equivalence transformations, which are in fact known as easily obtainable from 
elements of the corresponding equivalence algebra g~, can be factored out in the course of the 
computation. At the same time, the main ingredient of the automorphism-based version is 
finding the entire automorphism group of the algebra g~, which places a strong limitation on 
the dimension and/or structure of g~. This is why usually the automorphism-based version is 
applied to classes of differential equations whose equivalence algebras are low-dimensional. 

The megaideal-based version rests on the invariance of megaideals of equivalence algebras 
under the associated automorphisms and hence it does not require the explicit computation of 
the automorphism groups of equivalence algebras, which makes it suitable for finding the equiv¬ 
alence groups of classes with infinite-dimensional equivalence algebras. A disadvantage of the 
megaideal-based version is that the equivalence algebra g~ of a class of differential equations C\s 
may admit properties that induce essential constraints for equivalence transformations of 
but cannot be interpreted in terms of megaideals of g”^. 

Sometimes disadvantages of either of the versions of the algebraic method can be reduced by 
combining these versions. 
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We have restricted ourselves to the computation of the equivalence group of a single physi¬ 
cally relevant example, a class of nonlinear wave equations arising in nonlinear elasticity. The 
equivalence algebra of this class is infinite dimensional, which thus gives a prototypical example 
for the effectiveness of the megaideal-based method. We plan to write a more extensive paper 
with further examples on the use of both versions of the algebraic method in the future. 

Important further theoretical developments of the present method may include an extension 
to the computation of extended equivalence groups, which also play a role in the theory of group 
classification. The study of equivalence groupoids of classes of differential equations using the 
algebraic method will be another central milestone. 
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